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Abstract 

We study the angular power spectra of polarized Galactic synchrotron in the range 
10 ^ I < 800, at several frequencies between 0.4 and 2.7 GHz and at several Galactic 
latitudes up to near the North Galactic Pole. Electric- and magnetic-parity polar- 
ization spectra are found to have slopes around ole,b = 1-4 — 1.5 in the Parkes and 
Effelsberg Galactic-Plane surveys, but strong local fluctuations of (Xe,b are found 
at |6| ~ 10° from the 1.4 GHz Effelsberg survey. The Cpn spectrum, which is insen- 
sitive to the polarization direction, is somewhat steeper, being api = 1.6 — 1.8 for 
the same surveys. The low-resolution multifrequency survey of Brouw and Spoelstra 
(1976) shows some flattening of the spectra below 1 GHz, more intense for Ce,bi 
than for Cpn- In no case we find evidence for really steep spectra. The extrapolation 
to the cosmological window shows that at 90 GHz the detection of E-mode harmon- 
ics in the cosmic background radiation should not be disturbed by synchrotron, 
even around I ~ 10 for a reionization optical depth r r i <; 0.05. 
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1 Introduction 



Angular power spectra (henceforth, APS) of the Galactic synchrotron polar- 
ized radiation are raising up an increasing attention in these years. The first 
study, based on the Parkes 2.4 GHz survey (Duncan et al. 1997, hereafter 
D97), was due to Tucci et al. (2000), and other papers appeared shortly after 
(Baccigalupi et al., 2001a; Giardino et al. 2001). The main motivation for these 
works was the need for an angular-scale-dependent separation of the cosmic 
microwave background (CMB) signal from the polarized radio foreground; in 
fact the first work on the APS of a polarized Galactic foreground (Prunet et 
al. 1998) modelled thermal dust emission in view of the scenario to be met by 
Planck HFI at 143-217 GHz. CMB polarization is essential in order to remove 
degeneracies between important cosmological parameters (Zaldarriaga et al. 
1997), but the cosmological window (the region in the frequency-angular-scale 
plane where the cosmological signal is stronger than foregrounds) is narrower 
for polarization than for anisotropy. A careful study of Galactic contamination 
versus angular scale is thereby necessary. The study of synchrotron APS may 
also be important because of its bearing on the knowledge of Galactic struc- 
ture, and in particular of the transverse magnetic field in emission regions and 
the longitudinal field in compact foreground screens (Tucci et al. 2001). This 
point is open to future studies. 

The purpose of this work is to check the generality of the APS behaviour found 
by Tucci et al. (2000): Electric and magnetic parity APS are governed (with 
reasonable approximation for 100 < I < 800) by power laws Ce,bi l~ aE - B 
with oie,b — 1.4 — 1.5, in the portion of the Southern Galactic Plane probed by 
the D97 Parkes survey. Although such slopes are close to the values cc^^ ~ 
1.3 — 1.4 found for thermal dust emission (Prunet et al. 1998), later work 
cast serious doubts on the generality of this behaviour, as far as synchrotron 
is concerned. For the APS of the scalar PI = \fQ 2 + U' 2 somewhat steeper 
spectra, with api ranging from 1.7 ±0.2 to 1.9 ±0.3, are given for the Galactic 
Plane by Baccigalupi et al. (2001a) and Giardino et al. (2001). The former 
authors also find a much higher slope, ctp/ = 2.9 ± 0.2 for I ^ 10 2 , out 
of the Galactic Plane from the 1.4 GHz, low- resolution survey of Brouw and 
Spoelstra (1976: BS76). This discrepancy would be important when evaluating 
the synchrotron contamination at large scales; it would thereby be rather 
unconfortable in view of the role of CMB polarization harmonics with / ps 
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2 — 20 in the separation of effects of the primordial gravitational background 
from those of the secondary ionization of the cosmic medium. [This is one of 
the objectives of the SPOrt project, see e.g. Cortiglioni et al. (2001).] Note 
that the situation is not quite clear for the total intensity I either: Large 
values ai ~ 3 are supported by Bouchet and Gispert (1999) and Giardino 
et al. (2001), but oli ~ 2 follows from analysis of the Jodrell Bank 5-GHz 
interferometric survey by Giardino et al. (2000). 

Here we extend our analysis of APS in the Galactic Plane considering the 2.7 
GHz Effelsberg survey (Duncan et al. 1999: D99); further, we analyse several 
patches at intermediate Galactic latitudes, |6| < 20°, from the 1.4 GHz Effels- 
berg survey (Uyaniker et al. 1999: U99), and finally three regions from BS76, 
at 5 frequencies between 408 and 1411 MHz, covering latitudes up to near 
the North Galactic Pole. Following Tucci et al. (2001) and thereby extending 
the analysis of Baccigalupi et al. (2001a), we carefully distinguish between the 
APS providing a fuller statistical description of the spin-2 polarization field 
(i.e., Cei and Cbi, from which Cpi = Cei + Cbi is usually computed), and the 
Cpn spectrum which takes into account only the magnitude of the polarization 
pseudovector. This is necessary because Tucci et al. (2001) find quite signifi- 
cant differences between Cp\ and Cpn both for a Gaussian polarization field, 
as predicted in the standard scenario for CMB, and for synchrotron radiation. 

Our results confirm that slopes around 1.4 — 1.5 in the range 100 < I < 800 
are preferred on the average, in both the Northern and the Southern Galactic 
Plane for Cei, Cbi and Cpi, although significant fluctuations around theses val- 
ues are found in 10° x 10° patches. The fluctuations however become stronger 
in small patches out of the Galactic Plane, so that a preferred slope should 
not regarded as meaningful for the U99 survey. In fact, the partial regulari- 
ties that we find for polarization APS do not support the usefulness of global 
(i.e., full sky) Ce,bi for a satisfactory description of the spatial distribution 
of synchrotron: Local Ce,bi based on Fourier analysis are much more suitable 
to this purpose. A consistent picture, however, emerges at all of the scales we 
investigated. As a matter of facts, from the low-resolution BS76 survey we find 
moderate slopes for I < 70, at all frequencies and even near the Galactic Pole. 
Due to the limited range of I available for the fits, we find strong correlations 
between fitted parameters and, as a consequence, large error bars; however 
there is no evidence for steep spectra, since the best values lie in the range 
aE,B — 0.5 — 2.0 almost everywhere and are particularly small at frequencies 
where Faraday rotation is more important. For the Cpn spectrum, the slope 
fluctuations are smaller than for Ce,bi- Our best value for apj is 1.6 — 1.8 at all 
resolutions and Galactic latitudes, and all frequencies > 1.4 GHz. The differ- 
ence between Ce,bi and Cpn is significant, although not large in the average. 
The latter quantity should thereby be confronted with C P C ^f B \ not with the 
C , £ C £ / B ' ) spectra which are popular among cosmologists, and this is properly 
done in the last Section of this paper. 
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Table 1 



Sky regions used for computation of APS 



Ref. 


v (GHz) 


FWHM 


Galactic latitudes 


Galactic longitudes 


D97 


2.4 


10.4' 


-5° < b < 5° 


-122° < £ < 5° 


D99 


2.7 


5.1' 


-5° < b < 5° 


5° < £ < 74° 








4° < b < 20° 


45° < £ < 55° 








3.5° < b < 10° 


140° < £ < 153° 


U99 


1.4 


9.3' 


3.8° < b < 15° 


190° < I < 210° 








5° < b < 15° 


65° < £ < 95° 








-15° < b < -5° 


70° < £ < 100° 








-10° < b < 20° 


120° < £ < 180° 


BS76 


0.41 - 1.4 


2.3° - 0.6° 


26.5° < b < 62.5° 


4 = 60.5°, AV> = ±27° 








60.5° < b < 84.5° 


4 = 30°, Atj; = ±15° 



All of the results that we provide for polarization APS should not be affected 
by contamination from point sources; regions with flat Cu, on the other hand, 
are possibly dominated by point sources in total intensity, and are used in 
this paper to derive upper limits on their contribution to polarization APS, 
as discussed in Section 3. 



2 Data analysis 

2. 1 Polarized synchrotron surveys 

Our first high-resolution analysis of Galactic-Plane synchrotron polarization 
was performed on the 2.4 GHz Parkes survey (D97); this study, using twelve 
10° x 10° square patches (Tucci et al. 2000), is extended here to the 2.7 GHz 
Effelsberg survey (1999: D99) covering the additional Galactic longitude range 
5° < £ < 74° (see Table 1) and providing six more 10° x 10° patches. A belt 
extended about 54% of the Galactic Plane is thereby covered, although with 
two slightly different frequencies. The FWHM resolutions are 10.4' and 5.1' 
respectively, sufficient to achieve angular scales up to / ~ 10 3 . The nominal 
rms noise in D97 is 8 mK for total power and 5.3 mK for polarization (5.3 and 
2.9 mK in some more sensitive areas), while in D99 the rms noise is 9 mK. 

Our study is further extended to moderate Galactic latitudes analysing the 
maps of the Effelsberg 1.4 GHz survey (U99), consisting of five regions fairly 
close to the Galactic Plane, —15° < b < 20°. The angular resolution is there 
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Fig. 1. The selected patches in the BS76 survey. The patches are the same at all 
frequencies; the dots give the locations of the original BS76 measurements at 408 
MHz. 

9.35', and the rms noise is about 15 mK for total intensity and about 8 mK 
for linear polarization. The above regions having different sizes, we extracted 
different numbers of patches from them, namely, 1, 2, 2, 3 and 3 patches 
respectively, for regions corresponding to increasing row numbers in the U99 
sector of Table I. The patch size is 10° x 10° almost everywhere (so that a 
portion of the first region is not used), but only 6° x 6° in the second region. 

Finally we analysed the maps of Brouw and Spoelstra (1976: BS76) at five 
frequencies between 408 and 1411 MHz. The corresponding beamwidths and 
noise levels are summarized in Table 2. These maps cover a substantial (more 
than 40%) part of the sky, allowing us to compare results at quite different lat- 
itudes up to near the Galactic Pole, although with moderate resolutions. In or- 
der to get good sampling and large signal-to-noise ratios we selected three rect- 
angular regions, whose locations are reported in the last three rows of Table 1 
and in Fig. 1. As an example, we show the 820-MHz total polarization signal 
in Fig. 2. The BS76 maps are affected by undersampling, and the grid spacing 
(typically larger than the FWHM) depends on the sky region and frequency. 
In the selected patches the average grid spacing is about 2°. To make Fourier 
analysis feasible we thereby constructed three new, evenly spaced grids adopt- 
ing a Gaussian smoothing function with dispersion <r sm = FWHM/V8 In 2 
and FWHM = 3°. This procedure is straightforward near the Galactic Plane, 
where geodesies orthogonal to meridians can be identified with parallels with 
no appreciable errors. At high latitudes we constructed the grids projecting 
orthogonal geodesies, evenly spaced with Ab = 1.5° at starting points, out 
of the central meridian of each patch, having longitude £ = £ c as reported 
in the last two rows of Table 1. On such geodesies, extended respectively up 
to proper lengths = ±27° and ±15°, we picked out the centres of the 
Gaussian beams, with spacing Aip' = 1.5°. (Since these centres do not lie on 
parallels, the latitude ranges in the Table refer to central longitudes £ c .) We 
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Fig. 2. The total polarization (PI) field (in Kelvin degrees) for the BS76 survey at 
820 MHz, computed applying a Gaussian smoothing with F WHM = 3° to the Q & 
U fields. 
Table 2 



Parameters of the BS76 


survey 


v (MHz) 


FWHM (°) 


Noise (K) 


408 


2.3 


0.34 


465 


2.0 


0.33 


610 


1.5 


0.16 


820 


1.0 


0.11 


1411 


0.6 


0.06 



checked that for the chosen Aijj the proper distance between corresponding 
points on neighbour geodesies remains sufficiently close to Aip' even at the 
extremal longitudes i = £ c ± Aif). This in fact provides an intrinsic check for 
the applicability of flat-space concepts to celestial-sphere patches. 

2.2 Computation of angular spectra 

Because of the limited sky coverages, power spectra can be suitably obtained 
from Fourier analysis instead of the standard spherical-harmonic approach 
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(Seljak 1997). The estimators for the power spectra of fields can be derived 
by means of 



C 



XI 



(1) 



where Vt is the solid angle of the patch under analysis, iVj the number of Fourier 
modes in the interval around I chosen for averaging, w x l the noise contribution 
and bi the beam function. Note that X( I ) is the discrete Fourier transform of 
the smoothed field X{0i,<pi), a function of the bin coordinates which can be 
any of the Stokes parameters I, Q and U, but also (considering the electric 
and magnetic parity polarization) E and B. In the small scale limit, electric 
and magnetic modes can in fact be given in terms of Q, U modes by means of 
a simple rotation in /-space: 

E(l) = QCi) cos(20 t ) + 1/(7) sin(20 7 ), 

B(T) = -Q(l) sin(20 T ) + 17(7) cos(20 T ). (2) 

Thus dealing with maps of the fields Q and U, we end up with Cei and Cbi, 
which are most suitable for comparison with CMB spectra. Of course one could 
directly use Cq\ and Cjji, but these spectra are not rotationally invariant. The 
APS describing the total power of the polarization field (taking into account 
spatial variations of both magnitude and direction) is 

Cpi = Cqi + Cui = Cei + Cbi, I > 2. (3) 



Note the lower limit on /, arising from the spin-2 nature of linear polarization 
(namely, of the Q, U complex). Other authors have considered the APS of the 
polarized intensity PI = y/Q' 2 + U' 2 , which is defined for any I > 0. Equation 
(1) still applies with X = PI. The PI field however does not provide a 
complete description of polarization, being only related to the magnitude of 
the polarization pseudovector. The equality Cpn = Cpi is warranted only 
if the polarization angle is uniform inside all the survey area, and different 
spectral behaviours arise even for Gaussian random fields (Tucci et al. 2001). 

In conclusion as many as 7 APS can be defined, without counting cross- 
correlation and circular polarization spectra. Out of these, 4 are mutually 
independent, and out of the latter, 3 are polarization APS. 

The technique based on Eqs. (1) and (2), which has been already applied 
by Tucci et al. (2000) on D97 data, is fairly straightfoward for the small-size 
patches of high-resolution surveys. As done in that work, it is implemented 
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here with a cosine apodization to suppress border effects, and with subtraction 
of mean values to suppress aliasing. However, a careful analysis is required for 
BS76 maps. We have already remarked in the previous Subsection that for 
high-latitude patches we need geodesies to identify equally spaced grids to be 
used in Fourier analysis; now we observe that we must also take reference frame 
effects into account for Q and U fields. Such effects are especially important 
for our last patch under analysis, whose border is very close to the Galactic 
Pole. 

In particular, in order to smooth the Q and U fields we use the following 
procedure: (i) at each point within a proper distance x < 5<7 sm from a Gaussian 
beam centre, we construct the polarization vector P , (ii) we perform a parallel 
transport of P to the beam centre, obtaining 

pf T) = pg cos t + p<p sin t, pf T) = -p sin t + p<p cos t, (4) 



with t — J cos#d</>, and (iii) from p e ^ 'we recover the parallel-transported 
Stokes parameters Q( PT ) and U^ FT ^ and apply the weight function exp[— x 2 / (2crg m )] 
to them, obtaining the smoothed quantities 



(pt) = Efc Q ( k T) exp[-Xfc/ (2Qj 
^ sm E,exp[- X 2 fc /(2aL)] ' 

f/ (p T ) = Efc^ PT) exp[-^/(2q s 2 m )] 

E k eM~xl/(2^J] ■ { > 

We checked that in practice the parallel transport could be made with negligi- 
ble variations on any reasonable path close to geodesies, and a very convenient 
approximation is t — cos^ ave A0, with 9 me the average of polar angles at the 
starting and end points. The smoothed fields and £/ s ^ T ) are then used 

in Eq. (1) for the angular spectra of BS76 maps. 



3 Results 



3.1 High-resolution surveys 



Angular spectra for total intensity and E-parity polarization found in the 
Galactic Plane are shown in Fig. 3. For most of the square patches the curves 
are reasonably approximated by power laws, 

C xl = A x r°x. (6) 
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Fig. 3. The APS for total intensity (dotted lines) end E-mode polarization (solid 
lines), from 6 patches of the 2.7-GHz Effelsberg survey and 12 patches of the 2.4-GHz 
Parkes survey. 

This holds not only for X — I and E, but also for the B, Q, U, and P 
spectra which, having similar behaviours, are not shown in the Figure. In 
Table 3 we report the results of spectral fits for X = I, E and B, in the range 
100 <l < 800 for all of the 18 Galactic Plane patches. 

Total intensity spectra exhibit large variations in amplitude, by more than 
three orders of magnitude. (The variations appear magnified in the normal- 
ization parameter Aj, up to five orders of magnitude, because of a positive 
correlation with aj.) The same feature is not found in polarization spectra. 
This can be explained by observing that, while the total power emission de- 
creases fast when moving away from the Galactic centre, the polarized compo- 
nent remains much more uniform changing the Galactic longitude and latitude 
[see Figs. 4 and 6 in Duncan et al. (1995), and Fig. 9 in D97]; D97 noticed 
a "background component" in the polarized emission of about 20 mK, nearly 
constant over the entire survey. The distributions of the indices ax, too, high- 
light some differences between total intensity and polarization spectra: The 
values of aj range between 0.4 and 2.2, while the slopes of the polarization 
spectra in all patches remain relatively close to the mean value ~ 1.4 — 1.5. 
Moreover, low-emission regions show very flat spectra in total intensity, while 
no meaningful differences are found between high- and low-polarized emission 
regions. 
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Table 3 

Best-fit parameters for APS from D97 and D99 surveys 





Aj (K 2 ) 




aj 


A E (K 2 ) 


OLE 


a B 


250 


0.14x10" 


-3 


96 + - 12 
uyD -0.08 


0.13xl0" 4 


1.58 ±0.11 


1.45 ±0.11 


260 


0.20x10' 


-2 


1.48 ±0.12 


0.96xl0- 5 


1.48 ±0.12 


1.32 ±0.10 


270 


0.02 




1.02 ±0.11 


2.6xl0~ 3 


1.79 ±0.12 


1.25 ±0.12 


280 


0.48x10" 


-4 


u -++ -0.12 


0.44xl0- 4 


1.56 ±0.12 


1.78 ±0.11 


290 


1.8 




1.43 ±0.11 


0.32xl0" 3 


2.04 ±0.12 


1.90 ±0.12 


300 
310 
320 


0.07 

5.1 

0.48 




i 99 +0.12 
-0.08 

1.89 ±0.12 
1.76±0.11 


0.15xl0~ 4 
0.28 x 10~ 5 
0.40xl0~ 4 


1.56 ±0.12 
1.11 ±0.11 
1.49 ±0.11 


1.88 ±0.12 
96 + 012 

u - yD -0.09 

1.40 ±0.12 


330 
340 


0.21 
17.5 




1.46 ±0.09 
2.00 ±0.11 


0.22xl0~ 3 
0.48xl0~ 4 


1.57 ±0.12 
1.31 ±0.11 


1 74+0.15 

i -' + -0.09 

1.52 ±0.12 


350 


3.5 




1.64 ±0.11 


0.92 x 10~ 6 


u - 80 -0.06 


1.12 ±0.12 


360 


19.5 




1.65 ±0.12 


0.63xl0~ 4 


1.28 ±0.11 


1.33 ±0.12 


10 


3.40 




1.60 ±0.12 


0.21xl0~ 4 


0.92 ±0.10 


1 10 +0.08 
-0.12 


20 
30 


124.4 
11.5 




2.18 ±0.07 
1.79 ±0.10 


0.39xl0~ 3 
0.88xl0~ 3 


J-44_o.07 

1.50 ±0.10 


1 c-9+0.07 
L - oz -0.15 

1.63 ±0.10 


40 


0.10 




1.21 ±0.12 


0.25xl0~ 3 


1 36 + - 11 
L - 00 -0.06 


1 67+0- 05 
L - 01 -0.10 


50 


0.31 




i 90+0.21 


0.80xl0~ 4 


1.24 ±0.09 


1.52 ±0.13 


60 


0.65x10" 


-3 


1 00 + - 12 

i - uu -0.15 


0.71xl0~ 3 


1.68 ±0.07 


1.82 ±0.10 



^ Galactic longitude of the patch centre. The patch belongs to D97 for I = 
250° - 360°, and to D99 for I = 10° - 60°]. 

Table 4 



Average parameters from Galactic Plane fits 



Survey & method 


Ai (K 2 ) 


ai 


A E (K 2 ) 


OLE 


as 


D97, (a x ) 




1.37 ±0.44 




1.44 ±0.30 


1.46 ±0.29 


D97, (Q) 


2.2 


1.60 ±0.13 


0.12xl0~ 3 


1.53 ±0.11 


1.43 ±0.12 


D99, (a x ) 




1.71 ±0.43 




1.40 ±0.23 


1.57 ±0.19 


D99, (Q) 


10.3 


1.82 ±0.11 


0.31 xl0~ 3 


1.39 ±0.11 


1.55 ±0.12 


D97+D99, (ax) 




1.48 ±0.46 




1.42 ±0.30 


1.51 ±0.26 


D97+D99, (C e ) 


6.9 


1.72 ±0.10 


0.26xl0~ 3 


1.40 ±0.12 


1.54 ±0.10 



In Table 4 we report the values of the average slopes (rows labelled with 
(ax)), as well as the fit parameters from the average spectra (rows labelled 
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0.1 

lO" 2 
10" 3 
lO" 4 

O" 10 -6 

lO" 7 

10" 8 

10" 9 
0.1 

10~ 2 

10~ 3 

10- 4 
o~ 10 -° 

lO" 7 
10" 8 
10" 9 




b = [8,18] 
/=[45,55] 



b=[5,15] 
^=[65,95] 



b = [-5,-15] 
/=[70,100] 




100 



1000 



100 



1000 



100 



1000 



1 



1 



Fig. 4. The APS for intensity (dotted lines) end E'-mode polarization (solid lines), 
from 5 regions of the 1.4-GHz Effelsberg survey. 

by (Ci)) for both D97 and D99 surveys. The D97+D99 averages are computed 
rescaling the D99 data to 2.4 GHz using a spectral index (3 syn = —2.8 (Platania 
et al. 1998). The agreement is very satisfactory. On the other hand, the Cpn 
spectra turn out to be somewhat steeper; the average slopes turn out to be 
api ~ 1.6 — 1.8 (see Table 7 in the next Subsection), and fluctuations around 
these values are quite moderate. Although the difference with the above results 
is not large, it has a statistical significance, and is not surprising in the light of 
still large differences found in the arcmin angular range by Tucci et al. (2001). 



Fig. 4 and Table 5 report the results obtained from 5 intermediate-latitude 
regions of U99. (The E-mode spectra in the Figure again are illustrative of the 
behaviour of the other polarization modes.) The intensity spectra are found 
to be extremely flat (aj < 1) and low, indicating that the diffuse emission 
drops just out the Galactic plane. On the contrary, the polarization spectra 
do not show a decrease in amplitude with respect to the other two surveys; this 
means that the polarized background observed in D97 and D99 should extend 
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Table 5 

Best-fit parameters for APS from the U99 survey 



n°) 


H°) 


Aj (K 2 ) 




Otj 


A E (K 2 ) 




a E 


a B 


50 


13 


0.53x10" 


-5 


n QV+0-13 
u -' 5 ' -0.10 


0.16x10" 


-3 


1.22 ±0.08 


1.19 ±0.08 


143 


7 


0.35x10" 


-6 


+0.02 


0.22 




i n i a 
o k c; +0.14 

-0.17 


z - ' u -0.25 


150 


7 


0.50x10" 


-6 


Q +0.02 


0.80x10" 


-1 


oo +0.06 
z -' 50 -0.17 


9 in +0.23 
z - iu -0.16 


195 


10 


0.45x10" 


-6 


0. + - 02 


0.35x10" 


-1 


2.28 ±0.08 


2.32 ±0.08 


205 


10 


0.72x10" 


-6 


O.+o.n 


0.70x10" 


-2 


i n no 

i-99±o:?i 


1.98 ±0.10 


70 


10 


0.99x10" 


-4 


0.82 ±0.15 


0.73x10" 


-1 


2.41 ±0.11 


2.39 ±0.12 


80 


10 


0.89x10" 


-3 


1.13 ±0.13 


0.90x10" 


-3 


1.71 ±0.15 


1.79 ±0.19 


90 


10 


0.6xl0~ f 




o. + - 13 


0.11 




2.48 1°;^ 


2.23 ±0.13 


75 


-10 


0.16x10" 


-5 


0.16 ±0.13 


0.23x10" 


-5 


0.87 ±0.08 


1.50 ±0.12 


85 


-10 


0.12x10" 


-4 


0.87 ±0.12 


0.36x10" 


-5 


1.00 ±0.11 


1.17±0.12 


95 


-10 


0.20x10- 


-4 


0.81 ±0.11 


0.47x10" 


-4 


1.20 ±0.09 


fiS +0 - 07 

U- 0,3 -0.02 



at least up to £ ~ 15°. However, in U99 we observe large differences in the 
slope of polarization spectra from region to region: there are three patches with 
aE,B > 2 and two with q;_b,b ~ 1- Some of these regions cannot be considered 
as typical; for example, the area 140° < £ < 153°, 4° < b < 10.4° lies within 
the so called "fan region" , and the two regions centered at £ = 80° show rather 
complex structures. Such particularly large differences may be due to Faraday 
rotation effects, which are relevant at 1.4 GHz. However a significant dispersion 
in the results remains at higher frequencies, and undoubtedly shows how local 
(rather than global) spectra will be important for foreground subtraction in 
CMB polarization measurements. 

It is interesting to observe that several U99 patches exhibit nearly flat intensity 
spectra, i.e. aj ~ 0, with a remarkably uniform normalization. It can be 
argued that these patches are dominated (in total intensity) by extragalactic 
point sources, which are expected indeed to have flat APS as far as clustering 
can be neglected (Tegmark and Efstathiou 1996; Toffolatti et al. 1998). We 
can thereby estimate Cjf S ^ ~ 5 x 10~ 7 K 2 . To be quite conservative, we can 
take this number as an upper limit. From this limit, adopting a radio-source 
polarization degree of 5% [in agreement with De Zotti et al. (2000)] we get 
C^p < 1.3 x 10~ 9 K 2 for X = P, PI at 1.4 GHz. Assuming further a 
(frequency) spectral index /9rg — ~2 for radiogalaxies, we also get < 
1.5 x 10~ 10 K 2 at 2.4 GHz. We conclude that the contribution of point sources 
should be negligible in the whole range / < 800 for all of our polarization APS. 
This is consistent with our results: The reported Ce,bi, being ^ 5 x 10~ 9 K 2 for 
/ < 1000 at 1.4 GHz, do not show any flattening. Note also that more stringent 
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limits on (and therefore, on Cp'pu) might be derived from estimates of 

Tegmark et al. (1996) and Toffolatti et al. (1998). These however would be 
less safe. 

3.2 Low-resolution survey 




2 



log(l) 

Fig. 5. The APS for E-mode (solid lines) and B-mode (dotted lines) for the BS76 
patch centered at 6 = 5°. From top to bottom, the APS refer to frequencies of 408, 
465, 610, 820 and 1411 MHz. Curves labelled by a number n are shifted by a factor 
10 n . 

For the BS76 survey we analysed the maps at all of the frequencies, but due 
to the moderate resolution and the patch size, we could investigate only a 
limited range of I. In particular, uncertainties in the beam angular function 
for I > 7r/FWHM must exist in the original experiment; further, adopting a 
Gaussian shape for our smearing function exp[— x 2 /(2<7g m )] neglects the finite 
bin of the original measurements. Therefore, it is advisable to limit ourselves 
to the range such that b~[ 2 ^ e. In fact, the spectra computed by means of Eq. 
(1) turned out to increase for b^ 2 k, e, as can be expected for the difficulty of 
an accurate calculation of noise for the smoothed beams. Limiting ourselves 
to the range I < 70, this effect does not appear to visual inspection (see Figs. 
5-7); however, for the fits we adopt the modified function 

C xl = A x l~^+O x bi 2 . (7) 

Here O x , a parameter to be determined by the fit, is not intended to describe 
the field X, but rather to account for the inaccurate a priori estimate of 
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1 1.5 



logW 

Fig. 6. Same as Fig. 5 for the BS76 patch centered at b = 44.5°. 




1 1.2 1.4 1.6 1-8 2 



logW 

Fig. 7. Same as Fig. 5 for the BS76 patch centered at b = 72.5°. 

w x l . The upper limits on point source APS given in the previous Subsection, 
rescaled to the BS76 frequencies, exclude that the Oxb~ 2 term may mask 
important contributions from the above sources. As a matter of facts, for 
p RG = -2.3 we find for instance c£, S) < 4 x 1(T 7 K 2 for X = P, PI at 408 
MHz, some orders of magnitude below the measured APS. 

The results of the fits are given in Table 6. The declared uncertainties, which 
are 1-a errors on individual parameters computed from x 2 fields, are large be- 
cause of correlations in 3-parameter fits. In particular, we find strong, positive 
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correlations between Ax and ax, as illustrated by an example in Fig. 8, which 
refers to the 408-MHz low-latitude patch. This Figure shows x 2 iso-contours 
that we found in the (A x , ax) plane after minimization with respect to Ox- 
In spite of the above uncertainties, we find moderate slopes for a B and «b, 
generally in the range 1 — 2. At low frequencies, v < 610 MHz where Fara- 
day rotation is larger, we occasionally find some slopes < 1, and in one case 
(near the Galactic Pole in the 610 MHz map) a best value aE < 0. We should 
remark that due to the finite resolution of our sampling in (A x , ax, Ox) 
space, and the existence of narrow, elongated valleys with x 2 near minimum, 
the best values reported in Table 6 are less significant than the full extension 
of the confidence regions, described by the error bars in the Tables. However, 
none of the above considerations would change if the centres of such confi- 
dence regions were considered instead of the quoted minimum points. Thus 
there is sufficient evidence that polarization APS are somewhat flatter at low 
frequencies. 

The angular spectral behaviour does not exhibit any clear dependence on 
Galactic latitude at these moderate resolutions. It makes sense therefore to 
average over our three patches from the BS76 survey. Rows 1 to 5 of Table 7 
gives the average slopes ftp and apj for each BS76 frequency. At low frequen- 
cies, where Faraday rotation becomes more and more important, we observe 
a gradual flattening of both spectra; the deficit of steepness is more apparent 
for ap, in agreement with the considerations of Tucci et al. (2001). On the 
other hand, at 1411 GHz where Faraday rotation is less important (but not 
negligible at all), there is no evidence for any difference between ap and apj 
due to the large error bars. Also, the result quoted for apj at 1411 MHz is 
quite consistent to those found with higher resolution in the Galactic Plane 

10 1 P ' ' ' 1 ' ' ' 1 ' ' 1 ' ' ' E 




1 I i i i I i i i I i i i I i i i 

-20246 

Slope 

Fig. 8. Confidence levels for B-mode parameters fitted at 408 MHz for the BS76 
patch centered at b = 5°. The x 2 contours in the (Ap, ap) plane are reported at 
1-cr (solid lines) and 2-cr levels (dotted), after minimization with respect to Op. 
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Table 6 

Best-fit parameters for APS from the BS76 survey^ 



v (MHz) 


\ogA E 




log^B 


as 


log^p 


ap 


logA PI 


atpi 


408 


-0.9123 


i-ol2:S 


-o.7i;;| 


i 9+0.5 


-0 7+ - 1 
u - ' -0.1 


0.8l8i 


-o.2l8:i 


1 O + O.l 

i " 3 -0.3 


465 


o 1+0.8 


o.4l?:8 


-o.7l8i 


i r;+0.3 
l-°-0.9 


-0.9!?; 4 


1.212-5 


-0.518;? 


i 1+0.3 
1-1-0.2 


610 


-1.212;? 


i 9+0.4 
i - z -1.0 


_i n+ - 4 
1 - u -0.5 


i q+0.6 


-o.7l8:? 


i-i-8i 


-0.518;? 


l-°-0.3 


820 


i q+0.6 




i q+0.3 


1 4+ - 3 
l- 4 -0.3 


-o.9l8:? 


i 9+O.6 


-0.718;? 


i-5±8:I 


1411 


-2.oi;;| 


1 n+0.4 
i - y -0.2 


9 O+0.3 
_ ^-°-0.3 


i 9+0.4 
!- z -0.3 


-1.618;? 


1 S+0- 5 


-2.018;? 


i.8l8:5 


408 


-1-818:? 


i 9+0.4 
l- Z -0.7 


-1 5+ - 4 
i -°-0.8 


n 6+°- 6 

U.D_q 4 


-o.8l i:o - 


1-1-0.8 


_1 1+0-9 

i - L -o.5 


i q+l.l 
l-^-l.O 


465 


i n+0.6 
_J - z -0.5 


i r;+0.6 
i -°-0.8 


i 7+O.6 
— J -'-0.9 


i q+0.8 
i " 3 -0.9 


-1 0+0- 2 
l- u -0.8 


i 7+0.3 
J--' -1.1 


-1 0+ 1 - 


1 3+ 1 - 2 

!-°-0.3 


610 


-1 7+ 1 - 4 

1 -'-1.2 


1 3+ 1 - 5 


9 c;+0.5 
z -°-1.3 


1 l+°- 4 

i - i -0.9 


-i-9l2:l 


1 l+°- 4 


-1 2+ 1 - 6 


1 6+ L0 
l- D -0.9 


820 


9 q+0.5 
—2.1 


1 5 +1 -° 
l-°-1.8 


-2 0+ - 8 
z - u -1.8 


1 2+ L0 
l- z -1.8 


-i-8l?: 5 7 


i 1+0.9 
i - L -0.8 


_1 q+i-o 

i - y -1.3 


i 7 +i.o 

l-'-0.6 


1411 


9 c+0.3 
z - d -0.8 


1.4+°;? 


-1-418:5 


9 9+I.O 
z - z -0.8 


_1 ^+0-9 


1.91S:? 


-2.018;? 


i 7+0.9 


408 


-2.ol?:Z 


0.8l?;| 


-2.2123 


0.5±°;t 


_i 4+0- 4 

l- 4 -0.6 


0-418:! 


-0.4U;? 


i-il5:I 


465 


-l-Tlft 


n 6+ 1 - 4 

u -°-0.5 


-2.oi& 


o.5l8:l 


-1.712;? 


o.3l?;? 


-0.315:? 


o.9l°:I 


610 


O 7+1-2 

— J - ' -1.0 


-0 4 +L2 
u - 4 -0.8 


1 s+ 3 - 5 


i.ol8:8 


-1 4+ L4 


1 2+ 1 - 3 


-1 5+ 2 - 4 
i -°-0.1 


i 9+2.2 
l- z -0.9 


820 


9 1+0.8 
z -l-1.8 


1 2+°- 5 


_i 7+1-1 

i -'-2.1 


i c;+0.9 
l-°-1.5 


_i 4+°- 7 
i- 4 -o.8 


1 0+0- 5 


-1 5+ 1 - 2 
i -°-0.3 


i 9+O.8 
!- z -0.4 


1411 


_o n+0.5 
O.U_ 2 o 


1 4+ - 6 

1-+-2.1 


-1 3+ L0 
i " 5 -1.7 


9 i+l.O 


-1 5+ L4 
L -°-1.5 


1 9+ 1 - 1 
l- y -1.3 


-2 n+°- 5 

z - u -1.8 


i.8li:8 



1"Three sets of 5 frequencies refer to 3 patches, with increasing Galactic latitudes (cfr. Table 
1). 



Table 7 

Average slopes of Cp\ and Cpn spectra 



v (MHz) 


£-ran£ 


;e 


ap 


api 


408 


< 70 




0.56 ±0.24 


1.29 ±0.19 


465 


< 70 




1.05 ±0.43 


1.09 ±0.22 


610 


< 70 




1.10 ±0.14 


1.50 ±0.24 


820 


< 70 




1.14 ±0.28 


1.43 ±0.37 


1411 


< 70 




1.82 ±0.34 


1.76 ±0.51 


2417 


100- 


800 


1.48 ±0.12 


1.68 ±0.30 


2695 


100 - 


800 


1.50 ±0.11 


1.59 ±0.12 



at 2.4 and 2.7 GHz (cfr. the last two rows in the Table). 
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4 Conclusion 
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Fig. 9. CMB and synchrotron APS. The CMB Cei and Cpn angular spectra, com- 
puted for a sCDM model with a reionization optical depth r r ; = 0.1, are represented 
by the solid and short-dashed curves respectively. The corresponding synchrotron 
spectra are respectively given by the long-dashed and dotted lines; these are extrap- 
olated to 60 (upper lines) and 90 GHz (lower), assuming a frequency spectral index 
Asyn = —3. See the text for details. 

The regularities found for polarization APS in this work are not sufficient 
to establish the usefulness of global (i.e., full sky) C X i for a satisfactory de- 
scription of the spatial distribution of synchrotron itself, and a fortiori for the 
separation of Galactic synchrotron from CMB at higher frequencies. On the 
contrary, the local Cxi based on Fourier analysis are much more suitable for a 
fuller description of angular structure. This makes all-sky surveys of polarized 
synchrotron quite necessary. 

Quite significant fluctuations are found for parameters fitted in 10° x 10° 
patches. When we average over larger regions, the most stable behaviour is 
found for the polarized intensity spectrum Cpn which, for all of the surveys 
analysed by us at v > 1.4 GHz, everywhere exhibits a slope apj = 1.6 — 1.8. 
This holds in the full range 10 ^ I < 800, although we have large error bars 
for the range I < 70 investigated on BS76 maps. In spite of the more intricate 
situation for the other APS Cxi (X = E, B and P), we can state the following: 

• In the Galactic Plane, the slopes of electric and magnetic parity APS are 
quite moderate at 2.4 and 2.7 GHz, with averages ap ~ ap — 1.4 — 1.5 in 
the range 100 <l < 800. 

• Local fluctuations do not allow us to establish equally significant average 
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slopes out of the Galactic Plane in the same angular range. 
• At lower resolution, I < 70, large correlation between fit parameters cause 
large error bars; however, the best-fit slopes ap and ag stay in the range 
0.5 — 2 for almost all frequencies (in the range 0.4 — 1.4 GHz) and Galactic 
latitudes, and are quite inconsistent with values around 3. 

The above behaviours of polarization APS should be attributed to Galactic 
synchrotron, with no appreciable contamination from point sources. On the 
other hand, the total intensity APS may be locally dominated by sources, 
when they exhibit small amplitude and slope ai close to zero. 

Our results resolve the seeming discrepancy with other authors for the angular 
range 100 < I < 800, showing that investigators simply have to carefully 
consider which polarization APS is actually being computed. Our Galactic 
Plane result, a E — ap — 1-4 — 1.5, is very close to results found for the 
polarization APS of thermal dust (Prunet et al. 1998); this number is maybe 
deeply connected to Galactic structure. On the other hand, we do not confirm 
the high spectral slope found by Baccigalupi et al. (2001a) for Cpp at smaller 
I on three BS76 patches at 1.4 GHz. Our BS76 patches however are different 
from theirs, being expressely chosen to have a larger signal-to-noise ratio. Our 
best value apj = 1.76 ± 0.51, which arises from averaging over 1.4-GHz BS76 
patches, is consistent with our results at all resolutions and frequencies > 1.4 
GHz. The latter also agree with results obtained by Tucci et al. (2001) in 
the arcminute range. The slope of starligth APS, a^ 31 ' ~ 1.5 (Fosalba et al. 
2001), is closer to our a E ,B- 

Generally speaking, the PI and E & B fields contain different physical infor- 
mation. Since PI does not carry any information on the polarization angle, 
its spectrum cannot properly describe some related effects like, for example, 
beam depolarization; it should be used carefully, keeping in mind that it does 
not provide a complete description of the polarization field. 

In particular, if Cpp is extrapolated to the cosmological window, it is im- 
portant to make a proper comparison with the corresponding APS of CMB. 
Figure 9 compares both C E i and C Pn spectra of synchrotron and CMB. The 
reported CMB E-mode spectrum is the output of CMBFAST for a sCDM 
model with a reionization optical depth r r ; = 0.1. The corresponding Cpff B ^ 
spectrum is obtained through simulations in a 10° x 10° box, with the mean 
value being subtracted off. [Simulations of C^ MB ^ in the same box repro- 
duced the CMBFAST output to a satisfactory extent, see Tucci et al. (2000) 
for details.] The synchrotron APS in the Figure are extrapolated to 60 and 
90 GHz from the Galactic Plane average spectra (6th row in Table 4) for 
I > 100, and from BS76 patch centered at b = 44.5° for I < 70. In both 
cases, the computed normalization applies to high-emission regions and is not 
expected to be representative of the whole sky. For the extrapolations to high 
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frequencies we assume a spectral index f3 syIi = —3 (Platania et al., 1998). The 
results taken altogether offer a quite consistent picture. They imply that in 
high emission regions, whathever polarization APS is chosen, the synchroton 
polarized foreground should be comparable to CMB at 60 GHz even at large I. 
At 90 GHz the expected scenario looks more favourable for CMB experiments, 
both at small and large angular scales. Reionization effects on CMB should 
be investigated at I ~ 10 by means of C^ MB ^ rather than cj^f B ' ) . From the 
inspection of Fig. 9, and recalling that the height of the low-/ peak of C^ MB ' ) 
is approximately proportional r r 2 ;, we infer that the 90-GHz CMB signal from 
reionization should prevail on synchrotron at least for r ri <: 0.05. We finally 
remark that in the PI spectrum the "DC" signal C^ lB ^ might be interesting, 
but Fourier analysis is not appropriate to this purpose. 
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